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Abstract 

MPD thrusters have demonstrated between 2000 and 7000 seconds specific impulse at efficiencies 
approaching 40%, and have been operated continuously at power levels over 500 kW. These demonstrated 
capabilities, combined with the simplicity and robustness of the thruster, make them attractive candidates for 
application to both unmanned and manned orbit raising, lunar, and planetary missions. To date, however, only 
a limited number of thruster confiurations, propellants, and operating conditions have been studied. This work 
reviews the present status of MPD thruster research, including developments in the measured performance 
levels and electrode erosion rates, and theoretical studies of the thruster dynamics. Significant progress has 
been made in establishing empirical scaling laws, performance and lifetime limitations and in the development of 
numerical codes to simulate the flow field and electrode processes. 


Introduction 

Magnetoplasmadynamic (MPD) thrusters are currently being considered for applications as primary 
propulsion for both piloted and robotic and planetary missions. The high specific impulse provided by MPD 
thrusters, between 2000 and 7000 seconds, can reduce the spacecraft launch mass by more than a factor of two 
below that of chemical /aerobraked vehicles if the thruster systems operate with efficiencies over 50%. In 
addition, proper selection of mission scenario shows that the trip times can be made comparable to those 
achieved with high thrust systems (400 days to Mars) using reasonable extrapolations of current technology 1 * 3 . 
MPD thrusters have been tested in steady-state operation at power levels between 30 and 600 kW, and up to 20 
MW in pulsed mode for diagnostic purposes. The potential simplicity and robustness of the MPD thruster may 
be a significant advantage over alternative electric propulsion technologies. This paper presents an overview of 
the current status of MPD thruster technology, including both experimental and theoretical research efforts. 

Typical MPD thrusters, shown in Fig. la-c, consist of a central, thoriated tungsten cathode surrounded 
by an annular anode. The material used for the anode depends on whether it is actively or radiatively cooled. 
Propellant, usually argon or hydrogen, is injected into the rear of the chamber through an insulating wall. The 
propellant is heated, ionized, and accelerated by the current passing from anode to cathode and the Lorentz 
forces generated by the interaction of the current with either the self-induced or applied magnetic fields. The 
azimuthal self-induced magnetic field arises from the radial and axial distribution of the discharge current. An 
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external solenoid coil, typically placed outside the anode, is used to generate the axial and radial components of 
the applied magnetic field. The different electrode configurations, shown in Figure la - c reflect the diversity 
of the geometries studied: there are currently insufficient data to establish which thruster configurations will 
provide the required performance levels. 

A considerable amount of MPD thruster research and technology development was performed in the 
1960’s and very early 1970*s. Efforts were underway in government, industry and university laboratories to 
study both the self-induced and applied magnetic field devices. Between that time and 1988 most work focussed 
on purely self-field devices, and only in the last three years has there been renewed interest in applied-field 
thrusters. This renewal of interest in applied field thrusters arose primarily from the apparent inability of 
steady-state self-field devices to deliver the performance levels needed for propulsion applications. To date, the 
highest performance measured on a self-field thruster was obtained on a pulsed devise, which reached 5000 
seconds specific impulse at 40% efficiency with 1.5 MW input power 4 . By contrast, the highest performance 
applied field MPD thruster produced 7000 seconds specific impulse at 70% efficiency when operated in a 
steady-state mode on lithium propellant with 30 kW input power 5 . For propellants other than lithium, both self 
and applied-field steady-state thruster performance has been limited to efficiencies below 40% and specific 
impulses below 3000 seconds. In general, applied field thrusters have produced higher performance than the 
purely self-field devices. 

While the demonstrated performance levels indicate that MPD thrusters may be attractive for SEI class 
missions, a high power (200 kW - 1 MW), high performance, long life, MPD thruster has yet to be demonstrat- 
ed. MPD thruster life is currently limited by cathode erosion. The longest test to date, conducted in 1969, 
lasted only 500 hrs at 30 kW 6 . Efforts today are focussed on isolating the power loss mechanisms, identifying 
the dominant applied-field thruster acceleration processes, and minimizing electrode erosion. Direct testing of 
50 to 500 kW MPD thrusters for application to precursor missions 3 is underway, however, the stringent require- 
ments on the test facility preclude immediate testing of a steady-state multimegawatt MPD thruster in current 
space simulation chambers 7,8 . These facilities can be modified to increase pumping speed as described in Ref. 

7. The current facility limitations are forcing the development of scaling relationships and of experiments 
devoted to studying particular thruster phenomena, rather than a direct attack on high power thruster develop- 
ment. 


Two major reviews of the MPD thruster research and technology efforts have been performed. The 
first, by Nerheim and Kelly 9 , included both experimental and theoretical work conducted up to 1968, The 
second, by Sovey and Mantenieks 10 , focussed on experimental measurements of thruster performance and life up 
to 1987, though it did not include an assessment of the status of theoretical models. Over the past four years 
considerable progress has been made in improving the performance of 100 - 200 kW steady-state thrusters by 
using applied magnetic fields, and in identifying the performance loss and life limiting mechanisms for MPD 
thrusters operated at powers between 30 kW and 10 MW. In addition to an expansion of the experimental 
efforts, significant improvements have occured also in the numerical codes used to study MPD thrusters. These 
models have grown from one dimensional, inviscid, single fluid codes to fully two dimensional (axi symmetric) 
two fluid codes including viscosity, finite rate ionizatio,n and wall heat transport. In addition, several efforts 
are underway to determine the impact of plasma microturbulence on the plasma transport coefficients. 

The purpose of this report is to summarize the status of MPD thruster technology programs, including 
both experimental and theoretical components. The paper is divided into three sections, the first dealing with 
experimental measurements of thruster performance and power deposition, the second reviewing recent studies 
of MPD thruster cathode erosion, and the third presenting a summary of the ongoing theoretical work in the 
field. 
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MPD Thruster Performance 


The performance of any electric thruster is governed by the thrust efficiency, 17, and the specific 
impulse, 1^, which are calculated from: 


(a) 17 = 


T 2 

2mVJ 



( 1 ) 


where T is the thrust, m is the mass flow rate and V and J are the discharge voltage and current, repectively. 
Evaluating thruster performance therefore requires accurate measurements of thrust, propellant flow rate, and 
input power. By far the most challenging of these measurements is the thrust, which is subject to errors arising 
from the low thrust to weight ratio of the rocket, thermal and electromagnetic tares, and vacuum facility 
interactions. The recent work of Haag 11 illustrates the care that must be taken to ensure elimination of the 
thrust stand tares. Several previous efforts have addressed the impact of facility pressure on thruster perfor- 
mance 7 * 9 . Ambient pressure has been found to affect the thrust, terminal voltage, and electrode power 
deposition. These studies indicate that for open-throat thrusters (Fig. la) the test facility pressure during 
operation should be below 0.07 Pa, while for nozzle-type thrusters operating at high propellant flow rates this 
requirement can be relaxed to 1 Pa. In addition to thrust measurement errors, the actual mass flow rate in the 
thrust chamber may be significantly impacted by electrode erosion. This is particularly true for pulsed 
thrusters, in which the cathode mass loss may be 2-3 orders -of-magnitude higher than in a steady state thruster. 
In this work the data are reported as they appeared, no attempt was made to put them on a consistent basis with 
respect to potential sources of error. 

The highest values of efficiency and specific impulse reported over the last four years are shown in 
Fig. 2. The data were separated according to propellant, whether the thruster was operated with an applied 
magnetic field, and whether the thruster was operated in quasi-steady or steady-state mode. The numbers next 
to the symbols refer to the appropriate references 8,21,22,2630 . The power dissipated in the solenoid was not 
included in the performance calculations for applied-field thrusters. The highest performance reported, 40% 
efficiency at 3500 seconds specific impulse, was measured with a hydrogen propellant, quasi-steady, applied- 
field thruster. While no clear trends are evident, the only specific impulses over 2000 seconds were obtained 
using applied field thrusters. The performance measurements are replotted as a function of thruster power in 
Fig. 3 and 4. It appears from Fig. 3 that the efficiency does increase with input power, though there are 
substantial differences between the various thrusters tested. There is no such trend in the specific impulses 
shown in Fig. 4, as the highest value reported, 7000 seconds, was obtained using a steady-state applied-field 
thruster operated at 16 kW using hydrogen propellant. Note that the performance of the highest power steady- 
state thruster was rather low, reaching only 23 % efficiency at 950 seconds specific impulse. This low value of 
specific impulse is most likely a result of the high mass flow rate used for the test. For the operating conditions 
used, the thruster was in a predominantly thermal acceleration regime. 

While the applied-field MPD thrusters have demonstrated improved performance over self-field 
devices, little is known about their scaling characteristics. An empirical study where the anode and cathode 
radii and lengths were varied in a systematic fashion was recently performed to establish the required data base 1 . 
In that work, the anode and cathode radii were varied from 2.5 cm to 5.1 cm, and from 0.64 to 1.27 cm, 
respectively. The anode length was either 7.6 cm or 15.2 cm, and the cathode length was either 2 5 or 7.6 cm. 
A total of seven cylindrical geometries and one flared anode geometry were tested. The thrusters were all tested 
at an argon propellant flow tate of 0.1 g/s and a discharge current of 1000 A. The results of this study show 
that for a given mass flow rate and discharge current the thrust varies as 

r<x dbr 5+ * j/(I *> w 
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and the voltage as 


V <x 


Rl 

hRc 


B + k < f{L a ) 


( 3 ) 


where R, and R, are the anode and cathode radii, the k’s are constants, L ( is the cathode length, B is the 
magnetic field strength, and the second terms in these relationships simply reflect that while the anode length 
impacts both thrust and voltage, it does not affect the influence of the magnetic field. These results, and the 
data from which they were established, indicate that to increase the performance at a given applied-field 
strength, both electrodes should be short and that the anode to cathode radius ratio should be s mall However, 
it appears that while decreasing the anode radius increased the magnitude of the efficiency over the operating 
range studied, this may not remain the case for applied-field strengths over 0.1T. This is due to an observed 
increase in the rate of efficiency rise with increasing applied-field strength as the anode radius was made larger. 
While the impact of the propellant flow rate and discharge current were not quantified, comparisons of 
hydrogen and argon propellants show a substantial benefit from using the lighter gas. Figure 5 shows a 
performance comparison for the two gases, obtained with a cylindrical thruster having a 2.5 cm radius, 7.6 cm 
length anode operated on 25 mg/s of each gas at a discharge current of 750A (Fig. la). Not only does 
hydrogen provide approximately 1000 sec. higher specific impulse at a given field strength, but also the 
discharge is stable over a broader operating range, yielding almost a factor of two increase in ultimate 
performance. 

Performance Limiters 

The performance of MPD thrusters is limited by several factors, including electrode power dissipation, 
frozen flow losses, and the so-called "onset* phenomenon. The term "onset" is currently used to describe a 
wide range of phenomena, including high frequency voltage fluctuations, anode spot formation, and enhanced 
erosion of the chamber 12 ’ 13 . Over the past four years several efforts have been devoted to detailed studies of all 
of these performance limiters. 

The single largest power loss mechanism in MPD thrusters operated at power levels below 2-3 MW is 
anode power deposition 1 ' 14 " 14 . The power deposited into the anode is given by 

P a = J • + P C mv + Prad ( 4 ) 


where the first three terms arise from current conduction into the anode, P is a convective term due to the 
hot plasma flow, and accounts for radiation from the plasma. For low pressure devices (Fig. 1 a and c), 
the current conduction terms dominate, though for certain conditions P^, and P^, can reach 20% of the anode 
power 15 . For low power steady state thrusters up to 80% of the input power is deposited into the anode, leaving 
only 20% for plasma heating and acceleration. Several techniques for reducing the anode power fraction are 
currently under investigation, including injecting propellant near the anode surface and applying magnetic 
fields. The former technique increases the number of current carriers in the anode region and is thought to 
decrease V A> while adding a magnetic field increases the total power deposited into the thruster and decreases 
the fraction of that power deposited into the anode. The latter effect is illustrated in Fig. 6, which shows 
calorimetric measurements of anode power deposition for a 2.5 cm radius, 7.6 cm length anode cylindrical 
thruster. The data show that increasing the applied-field strength from 0.01T to 0.15T decreased the anode 
power fraction from 80 % to 60 % for the thruster geometry shown in Fig la. However, the total anode power 
increased linearly with increasing field strength. These results show that the magnetic field increases power 
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deposited into the plasma more rapidly than it deposits power into the anode. 

The cathode is not a significant power sink in steady-state thrusters, though its importance in quasi- 
steady thrusters, is poorly understood. Measurements made in steady-state thrusters show that the cathode 
consumes between 5% and 15% of the input power, depending on thruster geometry and operating 
condition 16 * 17 . No direct measurements exist for quasi-steady thrusters, though it appears that the cathode fall 
voltage (and thus cathode power deposition) is significantly higher for pulsed devices 18 * 19 . Terminal voltage 
measurements show a 20-40 V decrease when the thruster is operated steady-state as opposed to quasi-steady, 
with the most likely explanation being a decrease in cathode fall voltage 18 . This result not only indicates that 
power deposition measurements in quasi-steady thrusters may not apply to steady-state thrusters, but also implies 
that comparing theoretical voltage predictions with quasi-steady measurements may not be appropriate. 

The other MPD thruster power loss mechanisms, listed in approximate order of importance, include 
propellant ionization, heat transferred to the rear chamber wall, unrecovered plasma thermal energy, exhaust 
divergence, plasma radiation and atomic excitation. There is substantial evidence for thrusters operated on 
argon propellant that the plasma is nearly fully ionized for most operating conditions, leading to an ionization 
power sink of 3. 8x10* W/kg/s of flow rate. Ionization fractions for other propellants are not known. Prelimi- 
nary estimates of heat transfer to the rear chamber wall show that the dominant heat source is the radiation from 
the cathode surface, which is limited by view-factor considerations to less than one-half the cathode power. 
Estimates of the unrecovered thermal energy vary widely due to the large uncertainty in the ion temperature, 
which may be between 5000 and 100,000 K 20 . 

The ultimate limit on MPD thruster performance is a consequence of the long thruster lifetimes 
required for the missions of interest. The thrusters must be operated under conditions where the electrode and 
insulator erosion rate is low and where the risk of catastrophic failure due to arc instabilities is minimized. It is 
apparent from several studies 4, I2 * 20 ' 22 that there are certain critical operating conditions, hereafter called onset 
conditions", at which there is a rapid increase in the erosion rate and the appearence of localized, high current 
density arcs at the anode surface. The latter can lead to catastrophic failure of the thruster. For self-field 
thrusters both the onset condition and thruster performance are strongly correlated with the ratio P/ m, which 
fundamentally links the performance and life of these devices. It is interesting to note that several authors have 
found that using an applied magnetic field can extend the range of P/m over which a thruster can be 
operated 23 * 24 . In addition, applied-field thrusters may achieve the required performance levels at significantly 
lower discharge currents, permitting opration well below the critical values of P/m. A detailed review of 
empirical studies of the onset condition for self-field thrusters condition was completed recently by Preble 13 . He 
established correlations between the value of P/m at which onset occurred and the electrode lengths, the 
interelectrode spacing, propellant injection site, and propellant flow rate for cylindrical thrusters operated using 
argon propellant. The critical P/m increased when the shorter electrode (either anode or cathode) length was 
increased, and decreased when the interelectrode spacing was increased. Injecting propellant near the anode 
generally increased the critical P/m, though this trend also appeared to depend on the electrode geometry and 
magnitude of the propellant flow rate. Increasing the propellant flow rate decreased the critical P/m slightly. 
Preble was able to explain these correlations in a preliminary fashion by invoking an electrothermal instability , 
the threshold of which was reached when the propellant was fully ionized. 

Other experimental studies of the onset condition have focussed on plasma microturbulence. 

Researchers at Princeton University have experimentally measured the plasma dispersion relation and quantified 
the frequency dependence of density and temperature fluctuations in both 9 kW steady-state thruster and a 2-10 
MW quasi-steady thruster 30 * 25 . For both these devices, spanning three orders-of-magnitude in power level, they 
found that most unstable mode of plasma microturbolence was the lower hybrid drift instability. It is clear that 
the instability may have a dramatic and sudden effect on plasma transport properties, although the connection 
between the instability and thruster performance limits has not been quantified. 


5 



MPD Thruster Lifetime 


Cathode erosion has been identified as the major life limiter of MPD thrusters. First order calculations 
of propulsion requirements for lunar and Mars cargo vehicles indicate that the total impulse and life per thruster 
should be between 1x10 s to lxlO 9 N-s and 2000 - 8000 hrs, respectively. These estimates assumed a specific 
impulse of 2000-5000 s and a power range per thruster between 0.5 and 3 MW. The lowest, most reliable 
cathode erosion rates measured prior to 1987 yielded predicted total impulse values that were two to three 
crders-of-magnitude below the estimated requirements defined above 10 . Such estimates, however, were difficult 
to make because the reported erosion rates varied over five orders-of-magnitude. Within the last few years 
significant progress has been made in both the measurement of erosion rates, cathode surface diagnostics, and in 
the theoretical understanding of MPD thruster cathode erosion processes. 

The most recent measurements of steady-state and quasi-steady thrusters are listed in Table I, where the 
erosion rates are listed as mass lost per unit charge transferred through the cathode surface. The mass loss rates 
range from 5X10* 4 to 0.2 /xg/C for steady-state thrusters and 0.2 to 60 pgfC for quasi-steady thrusters. Clearly 
quasi-steady thrusters have much higher cathode erosion rates than steady-state thrusters, though there is some 
overlap for quasi-steady thrusters operated with low work function, barium impregnated, cathodes. The lowest 
values of steady-state thrusters were obtained at the University of Stuttgart using a nozzle-type thruster operating 
with high purity propellant. 

The large differences between quasi-steady and steady-state erosion rates are readily explained by 
considering the electron emission mechanisms for the two operating modes. For quasi-steady thrusters, electron 
emission from a cold cathode is achieved by field enhanced thermal emission. The current attaches the cathode 
via highly mobile arc spots, typically between 10 and 50 pm in diameter 31,36,38 , characterized by high current 
density (> 10 11 A/cm 2 ) and high local temperature (> 5000 K). This emission mode is characterized by 
tungsten droplet ejection and strong evaporation, which is manifested by small "craters" from several microns to 
a tenth of a millimeter in size 31 . The cold cathode erosion rate of quasi-steady thrusters depends on the pulse 
frequency, pulse duration and the power applied. The use of an applied magnetic field reduced the cathode 
erosion rates by approximately a factor of two in two of the tests reported in Table l 23,29 . 

Several investigators have varied the duration of steady-state tests to permit separation of the erosion 
rates characterizing the ignition and steady-state phases of operation. The ignition mode erosion rate generally 
varies from 5-30 ^g/C 1718 * 33,39 , values which are close to those measured for quasi-steady thrusters operated with 
the same cathode material. The lowest erosion rate during the ignition mode was measured during single 
"shot" measurements. Higher erosion rates were measured during multiple "shot" runs, with the erosion rate 
decreasing with increasing number of "shots" during a test 33 . The type of propellant did not effect the ignition 
mode erosion rate 33 . For steady-state MPD thruster operation lasting many hundreds of hours the ignition rate 
should not significantly contribute to the total erosion rate of the cathode. 


The steady-state thermal emission mode is reached when the cathode reaches a sufficiently high bulk 
temperature to permit electron emission in a diffuse mode of current attachment. Surface temperature 
measurements made using imaging techniques 17 show that thermionic emission governed by the Richardson 
equation appears to explain electron emission from the cathode 17,33,35 . This transition from ignition mode to 
steady state diffuse mode has been observed using high speed photography 17 . During the steady-state mode the 
surface temperature of the cathode can be substantially below 3200 K with current densities between 10 - 100 
A/cm 2 . 


The relatively large differences in the steady-state erosion rates listed in Table 1 can be explained by 
considering the erosion mechanisms and the different thruster geometries used for the tests. Simultaneous 
measurements of cathode erosion rates and surface temperature distributions show that evaporation plays a 
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dominant role in the cathode mass loss. Estimates show that sputtering erosion does not contribute significantly 
for typical cathode temperatures and ion energies. While the maximum erosion rate appears to be determined 
by the vacuum evaporation rate of the cathode material, the actual net mass loss rate can be dramatically 
reduced by backscattering of the evaporated atoms off ambient gas atoms which return the material to the 
cathode surface. In addition, ionizing collisions between evaporated atoms and electrons yield ions which are 
pulled back to the cathode by the electric field. The influence of these processes on the cathode mass loss rate 
can be estimated by evaluating the appropriate mean-free paths as suggested by Myers 17 , or by considering the 
ambient ^apor pressure as suggested by Polk 35 . The hypothesis that the erosion rate of a cathode is governed by 
the backscattering mechanisms was confirmed at Stuttgart 33 and in the Soviet Union 37 . As shown in Fig. 7 for a 
bench top arc experiment, the erosion rate was found to decrease with increasing ambient pressure. Indirect 
evidence of this phenomena is the constant erosion rate measured as a function of discharge current by Schrade 31 
at low discharge current levels (2-4 KA) and the relatively large variation of cathode erosion rates measured as 
a function of discharge current by Polk 36 . It is clear that the influence of ambient pressures on backscattering 
rates fundamentally couple the thruster life, the thruster geometry, and the operating conditions. 

For the nozzle-type thrusters operated at the University of Stuttgart, the ambient pressures sorrounding 
the cathode lead to high backscatting rates. In this case, the net loss due to evaporation is reduced, and other 
mass loss mechanisms play a larger role. This may explain the impact of propellant purity found at the 
University of Stuttgart, which indicates that surface oxidation plays a significant role 18,33,35 . These tests, 
conducted with hydrogen, argon, and nitrogen, clearly showed that eliminating water and oxygen from the 
propellant using a chemical filter reduced the cathode erosion rate by up to a factor of twenty. This is the first 
time this effect has been reported for MPD thrusters. 

i 

The importance of the cathode material arises due to its influence on the surface work function, which 
is the energy required to liberate an electron from the metal. A small reduction in work function can lead to a 
large drop in the cathode surface temperature required to emit a given current density. This, is turn, can lead 
to a large drop in the material evaporation rate as it is exponentially dependent on temperature. For this reason 
most cathodes utilize a small amount of low work function impregnant material, usually 2% thorium oxide. 
Experiments have clearly shown the importance of having the low work function material 33 . However, for the 
thoriated tungsten the evaporation rate of the thorium from the surface exceeds the diffusion rate to the surface 
when the cathode temperature exceeds 1800 to 2300 K 17,35 . Since the surface temperature of a 2% thoriated 
tungsten cathode in an MPD thruster usually exceeds these limits, the net loss of the low work function material 
will result in a relatively quick failure of the cathode, unless backscattering limits the escape of low work 
function material. Experiments using both new and slightly used (—2 hrs) cathodes indicated a small Th0 2 loss 
during the first several hours of operation 33,35 . Several hours of operation after the "conditioning 1 ’ of the 
cathode indicated no more Th0 2 loss as inferred from a constant tip temperature during the test 35 . Some recent 
tests have explored using barium oxide impregnated cathodes in steady-state thrusters 32 though the results were 
inconclusive. 

The discussion of the erosion mechanisms and the results of recent erosion rate measurements of an 
MPD cathode operating in a steady-state diffuse mode leads to the following conclusions: 1) Evaporation 
appears to be the dominant erosion mechanism of MPD thruster cathodes, 2) The ambient pressure at the 
cathode surface can substantially affect the erosion rate of the cathode, 3) the impact of propellant purity 
appears to be substantial for thrusters operated at relatively high pressure, 4) the loss of low work function 
material impacts the cathode lifetime via its influence on the surface temperature. These factors indicate that the 
simple loss parameter of mass loss per unit charge is not adequate to evaluate erosion rates of MPD thruster 
cathodes. 


MPD Thruster Models 

As discussed in the previous sections, MPD thrusters have undergone extensive experimental 
development since the early 1960’s. In contrast, a detailed theoretical understanding of MPD thruster operation 
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has been slowed by the complex nature of the coupled electromagnetic and gasdynamic acceleration processes. 
Simplified analytical and numerical simulations using ID and quasi-ID approximations of the fluid 
magnetohydrodynamic equations have provided valuable insights into both self-field and applied-field MPD 
thruster operation, but they are constrained in their ability to predict detailed thruster performance. The 
removal of such limitations via 2D and quasi-2D modeling has recently become practical with the emergence of 
high-speed computational facilities, allowing model validation and refinement using the existing experimental 
data base, while in return establishing a theoretical basis to guide further experimentation. 

Governing Equations 

The equations used to model MPD thruster dymamics are generally based on the fluid 
magnetohydrodynamic (MHD) equations, which are valid in the collisional plasma regimes of interest for MPD 
thruster operation 40 . The full set of equations includes coupled electromagnetic field and plasma flow equations, 
which together with expressions for the equation of state, electromagnetic and thermal transport coefficients, and 
ionization/recombmation processes provide a complete description of the system. A fundamental exposition of 
the MHD formulation is provided by Sutton and Sherman 41 , the salient features of which are reproduced below. 

Electromagnetic Field Equations. 

The basic set of electromagnetic equations includes the full complement of Maxwell’s equations: 

(a) V • B = 0 (6) 

(c) V • ^ = pj( 0 ss o (d) 


V x B = no J 

V x E = -~ 

at 


(5) 


where E and B are the vector electric and magnetic fields, respectively, and j is the electric current density. 
The current density is related to the electromagnetic fields and plasma flow properties through a generalized 
Ohm’s law, usually given in the form: 

J = . [i + (? x B']] - ^ (J x i) + (1 - a) J ^ [(J x I) xb ] (6) 


The quantity a is the electrical conductivity, v is the plasma flow velocity, and a is the ionization 
fraction of the plasma, fi is the electron Hall parameter, which is the product of the electron cyclotron 
frequency and the electron collision time, while Oj is an equivalent expression for the ion Hall parameter, 
defined in terms of the ion cyclotron frequency and ion-neutral collision time. The ion Hall parameter enters 
into the ion slip term, which allows for incomplete coupling between the ions and neutral gas particles in a 
partially ionized plasma. 

The electric field distribution may be obtained by rearranging Ohm’s law: 


E = -V<f> = - 
a 


- n /- - 
1+b(i* b 


)]-(»xB) 


( 7 ) 


where ^ is the electric potential. The radial component of the electric field may then be integrated from 
cathode to anode to determine the potential drop across the plasma. 
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A magnetic transport equation may be derived by combining Maxwell's equations with the generalized 
Ohm's Uw: 


V x E = 



(< V » *)] 



( 8 ) 


Solving the azimuthal component of the above equation provides the induced azimuthal magnetic field 
distribution. The azimuthal magnetic field is a function of the total discharge current J, hence a streamfunction 
r(J/r) — J may be used to represent lines of total enclosed current. 


The radial and axial current densities may be obtained from the calculated self-magnetic field 
distributions: 

J r = -■ 


1 BB e 


Mo dz 

* " iilrrW) 


( 9 ) 


The azimuthal current density j # , which is due to the applied magnetic field, may be obtained from the 
azimuthal component of the general Ohm's law. The applied magnetic field components may be specified a 
priori or obtained from the solution of Maxwell's equations for an external current distribution. 

Fluid Equations 

The fluid equations encompass the conservation equations of mass, momentum, and energy, coupled 
with thermodynamic and caloric equations of state, ionization and recombination processes, che mical kinetics if 
the flow is reacting, and expressions for the plasma transport coefficients. 

The conservation of mass for a compressible fluid is given by: 


^ -f V • OnT) = 0 (10) 

where p is the plasma mass density. Excluding ambient gas entrainment or electrode erosion, the mass 
conservation law states that the amount of propellant mass injected into the thruster per unit time must equal the 
amount of mass exhausted from the thruster per unit t ime . 

The conservation of momentum yields an equation of motion for the plasma, given by: 

'’(£0 = ~ 7 '’ + ( ,,x ®) + ? ( ii ) 

where p is the plasma pressure, if is the viscous body force, and D/Dt denotes the advective derivative: 


D_ 

Dt 


! + (»-• 7) 


( 12 ) 
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For a self-field thruster, the components of the electromagnetic j xB body force simplify into an axial 
(J,B 8 ) body force, which provides direct electromagnetic thrust ("blowing"), and a radial (-j 4 B 0 ) body force, 
which provides electromagnetic compression of the plasma and a subsequent pressure force along the cathode 
surface ("pumping"). The electromagnetic body forces become significantly more complex for an applied-field 
thruster, where the interaction of the discharge current with the applied magnetic fields give rise to azimuthal 
body forces. The azimuthal body force can impart rotational motion to the plasma, which might be partially 
recovered as thermal energy (leading to enhanced electrothermal thrust) or converted to directed thrust via 
interactions with the applied magnetic field. Additional radial and axial body forces may be generated through 
interactions between the induced azimuthal current density and the applied magnetic field components. Viscous 
body forces and unbalanced pressure gradient forces must also be incorporated into a complete description of 
the plasma flow. While experimentally robust, the MPD thruster operates via some rather complex acceleration 
mechanisms. 

A general form for the conservation of energy, which expresses a balance between the energy sources 
and sinks in the plasma, is given by: 

p Wt = m + + * x S) ' v ‘ f+ * v (13) 

The quantity h is the specific enthalpy, where the second term on the right is the Joule heating term, 
is the viscous dissipation function, and q is the heat-flux vector. For models' incorporating more than one fluid 
species (e.g. ions, electrons, and/or neutrals treated separately), individual conservation equations for each 
species must be used. The species momentum and temperature equations are then coupled through collision 
integrals to allow interactions between each plasma component. Ionization and recombination processes may be 
used to evaluate the ion, electron, and neutral species number densities, in concert with a quasi-neutrality 
approximation such that the calculated ion and electron number densities are equivalent. 

The fluid equations are closed with a thermodynamic equation of state relating the plasma pressure, 
density, and temperature, a caloric equation of state relating the specific enthalpy to the plasma density and 
temperature, and expressions for the plasma transport coefficients (viscosity, thermal conductivity, and electrical 
conductivity). 

Approximations 

The full set of coupled, nonlinear partial differential equations used to model the MPD thruster plasma 
interactions present a formidable challenge for current numerical simulation techniques. Consequently, authors 
have employed simplifying assumptions to make the system dynamics more tractable. Table 2 lists the 
approximations used in several recent MPD thruster models 42 ^ 53 , defined according to the assumed plasma 
properties and model geometries. The utility of the various approximations, and their validity for modeling the 
MPD thruster, are discussed below. 

Plasma Species 

The fluid approximations outlined above are valid for a coliisional plasma, in which the electrons and 
ions separately relax to local thermodynamic equilibria over time scales which are short compared with changes 
in the overall plasma conditions, and in regions small compared with the overall plasma dimensions. The most 
simple fluid approximation is to ignore differences in the electron and ion species, and assume the plasma 
consists of a single conducting fluid which is described by a single mass density and a single (bulk) fluid 
velocity. This one-fluid. approach leads to some ambiguity in definitions of the plasma transport coefficients, 
which are usually defined in terms of the separate electron and ion transport properties 41 . In addition, such 
effects as ambipolar diffusion, electrode sheath formation, ionization processes, and plasma instabilities cannot 
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be directly modeled with a single fluid approximation. Using a multi-species model (separate ion and electron 
species, and possibly an additional neutral species) circumvents these limitations, and provides a more accurate 
treatment of the transport phenomena; however, additional accuracy is obtained at the cost of increased 
complexity in solving separate density, momentum, and energy equations together with the associated species 
momentum and energy coupling models. Due to this added complexity, the one-fluid approach has been used 
extensively in MPD thruster modeling, and has provided a somewhat crude but effective method to gain insight 
into the general plasma behavior under a variety of operating conditions. 

Plasma Temperature 

Several of the models listed in Table 2 have employed a single temperature approximation in 
conjunction with a single fluid assumption to simplify the governing equations. While providing a useful first 
estimate of the plasma dynamics, the single temperature approximation may lead to gross inaccuracies in the 
plasma transport coefficients, which are sensitive to the plasma electron and ion temperatures 41 . Models which 
incorporate separate species energy equations (Table 2) have predicted significantly different electron and heavy 
particle temperatures. Such results indicate the plasma may not be in thermal equilibrium in many operating 
regimes of interest, and illustrate the disparity between single fluid models and the real plasmas which they 
model. Although the complexity of a multi-species, multi-temperature plasma model may not be warranted in 
the initial stages of model development, it is clear that separate electron and ion temperatures should be 
incorporated to provide more accurate predictions of the plasma transport coefficients. 

Plasma Ionization 

The assumption that the plasma is fully, singly ionized has been incorporated into most single fluid and 
multi-species MPD thruster models, despite the varying degrees of ionization observed experimentally 16,64 . This 
assumption does not significantly affect the classical plasma transport coefficients due to their weak logarithmic 
dependence upon the plasma number density. However, the transport coefficients are sensitive to the plasma 
species collision times, which in turn are inversely proportional to the plasma species number desities. 
Consequently, an accurate calculation of the plasma ionization fraction is required for an accurate evaluation of 
the Hall terms appearing in the generalized Ohm's law. In addition, a knowledge of the ionization state of the 
plasma is important for accurately modeling electrode sheaths 65,66 , boundary layer flows 67 , and plasma 
instabilities I3 * 30 - 6B - 75 , which are not generally included in MPD thruster flow models due to the added complexity 
they bring to the problem. Discussions of nonequilibrium plasma ionization processes in self-field MPD 
thrusters may be found in references 76-80. 

Miscellaneous Plasma Properties. 

Hall Effect . 

In the absence of a magnetic field, conduction electrons will follow electric field lines between the 
cathode and anode. In the presence of either induced or applied magnetic fields, these simple electron 
trajectories become curved 41 , giving rise to a transverse current ("Hall" current) in a direction orthogonal to 
both the electric and magnetic fields. The magnitude of the transverse Hall parameter (0) depends upon the 
ratio of the electron cyclotron frequency to the electron collision frequency. The Hall current becomes 
important when fl * 1 , and dominates the system dynamics when 0 > > 1 . The Hall effect skews the current 
distributions between the electrodes, causing current to concentrate upstream along the cathode and downstream 
along the anode (Fig. 8). 

Viscosity. 

Viscous effects cause a dissipation of directed kinetic energy into thermal energy with a consequent 
reduction of thrust efficiency. Viscous dissipation may significantly contribute to heating the heavy plasma 
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species, driving the ion temperature to equal, and in some instances exceed, the electron temperature 5 *** 0,11 . 
Jon Slip. 


All single fluid and most multi-species models assume the ion motion is equivalent to the bulk plasma 
motion. In operating regimes characterized by low plasma densities and high magnetic field strengths, however, 
the collisions between plasma ions and neutral particles may be insufficient to keep the ion motion coupled to 
the neutral gas flow. The ions may then cross the interelectrode gap or proceed on their own tiajectories 
without contributing to the bulk fluid acceleration, causing a decrease in thrust efficiency. In single fluid 
models, and in models which assume a fully ionized plasma, the ion motion is identical with the bulk plasma 
motion, and the ion slip term is zero by definition. 

Status of MPD Thruster Modeling 


In their review paper, Nerheim and Kelley* provided an account of the general status of MPD thruster 
theory as of 1967. Relevant sections of their review article are summarized below, to provide a benchmark 
against which more recent MPD thruster modeling efforts may be compared. 


The electromagnetic force produced by the interaction of an applied current and self-induced magnetic 
field had been derived analytically by Maeckei* 2 in 1955, and incorporated to explain self-field MPD thruster 
acceleration processes by the mid- 1960’s. The electromagnetic force, in the absence of an applied magnetic 
field, has the form: 



(14) 


where F, is the self-induced electromagnetic force (thrust), J is the total discharge current, r t and r c are the 
anode and cathode radii, respectively, and C is a constant, which ranges in value from 0.25 if the current 
distribution is uniform over the cathode surface (excluding the cathode tip), to 0.75 if the current is uniformly 
distributed over the cathode tip* 3 . This expression is popularly known as the "bJ 2 law", where b incorporates 
the various constants into a single thrust coefficient. The bP law works surprisingly well to predict self-field 
thrust when the thruster geometry is sufficiently simple (such that r, and r c are well defined) and electrothermal 
thrust contributions can be ignored. The simple derivation fails for thruster geometries more complex than 
cylindrical electrodes, and for thruster operating regimes where the thrust is produced primarily by thermal 
expansion of the propellant rather than electromagnetic acceleration. 

Although thrust-producing mechanisms have been postulated for applied-field MPD thrusters, no 
comparable analytic expression relating electromagnetic thrust to engine operating parameters has yet been 
derived. One such mechanism discussed in Nerheim and Kelley (and the references therein) invokes the 
electromagnetic interaction between the applied current and the applied magnetic field to create an azimuthal 
torque on the plasma, which imparts rotational kinetic energy to the ionized propellant. The rotational energy 
might then be converted to directed kinetic energy through expansion in a magnetic nozzle* 4,15 , since the angular 
velocity must decrease to conserve angular momentum as the plasma expands, and the decrease in angular 
velocity must be accompanied by an increase in axial velocity to conserve kinetic energy. A second applied 
field acceleration mechanism requires the existence of Hall currents* 6 . Applied magnetic fields might induce an 
azimuthal Hall current in the plasma, which in turn interacts with the radial and axial components of the applied 
magnetic field to produce axial (-j,B f ) and radial (J^B,) body forces, which contribute to directed thrust and 
radial compression of the plasma. The experimental evidence for the existence of Hall currents in applied-field 
MPD thrusters was limited and contradictory, however, and their role in applied-field MPD thruster acceleration 
was uncertain*. 
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Summarizing the status of MPD thruster theory in 1967, Nerheim and Kelley concluded that the level 
of understanding of basic MPD thruster operation was insufficient to "...allow a critical evaluation of the 
reported performance data, to extrapolate the performance obtained in the laboratory to that which may be 
expected in space, or to aid appreciably in the design of an ’optimum’ engine. Thus, the fiill potential and the 
limitations of the MPD thrustor [sic] are not known at the present time". While the latter comment remains 
true, steady progress has been made in several areas of MPD thruster modeling and analysis. The status of 
recent MPD thruster modeling efforts, and associated issues which must be resolved, are outlined in the 
following sections. 


Self-field MPD Thruster Models. 

Several of the self-field models listed in Table 2 have been designed to investigate performance trends 
rather than provide specific calculations for comparison with experimental measurements. While limited in 
scope and detail, these models have provided necessary insight into several aspects of self-field MPD thruster 
operation. A few of the models have made specific predictions of such global parameters as thrust, current 
distribution, and voltage, which can be compared with experimental measurements. These specific comparisons 
and more general trends are discussed below. 

Thrust. 


Several recent MPD thruster models have been fairly successful in predicting self-field MPD engine 
thrust and thrust coefficients over a wide range of operating conditions 47,55 * 57,62 . Figure 9 presents measured 87 
and predicted 62 thrust characteristics for the Princeton University half-scale flared anode thruster (HSFAT), 
plotted as a function of discharge current for various argon propellant mass flow rates. Good agreement has 
been obtained between the NASA LeRC single-fluid, two-dimensional model and the experimentally measured 
values for each of the mass flow rates and discharge currents up to F/m values of approximately 1.5 x 10 n A 2 - 
s/kg, at which point the calculated thrust begins to deviate from the measured values. Thruster voltage 
oscillations were experimentally observed at F/m values around 1.8 x 10 11 A 2 -s/kg, and it is probable that the 
simple single-fluid model used in the comparison cannot accurately model the physical processes occurring in 
the thruster as this operational instability is approached. Good agreement between predicted and measured 
thrust has also been obtained by the University of Stuttgart 55 for the DT-series thrusters operated at low values 
of F/m. 


Figure 10 reproduces experimentally determined values of the thrust coefficient for the HSFAT as a 
function of F/m with predictions from the LeRC model overlayed on the experimental data 62 . Again, good 
agreement is observed between the model and the experimental measurements. Similar agreement has also been 
obtained using an analytic self-field model developed by Martinez-Sanchez 47 . The thrust coefficient varies from 
a value of approximately 7.5 x 10 7 N/A 2 at low F/m, corresponding to a predominantly electrothermal regime 
of operation, to an approximately constant value of 1.1 x 10* 7 N/A 2 at values of F/m which correspond to 
predominantly electromagnetic regime of thruster operation. 

Current Distribution . 

Predictions of steady- state anode current distributions have compared favorably with measurements 
performed at the University of Stuttgart using cylindrical, segmented anodes 51,55,62 . Table 3 displays the 
experimentally measured anode current fractions for the ZT-1 MPD thrustei* 8 , with the current fractions 
predicted by the Stuttgart 51 and LeRC 62 2-D MPD models. Good agreement is obtained between the models and 
the measured anode current fractions for the operating data reported (6000 amps, 6 g/s of argon). The 
predicted current distributions within the thruster are displayed in Fig. 11, showing similar distributions between 
the two models. Current distributions within the steady-state thruster could not be measured, but the agreement 
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with the anode current fractions lends some confidence that the predicted steady-state distributions are 
reasonable. 


This confidence quickly erodes, however, when predicted contours are compared with actual current 
contours measured in quasi-steady MPD thrusters. Figure 12 displays predicted 62 and measured* 7 current 
contours for the Princeton HSFAT, operated at 17.8 kA with an argon mass flow rate of 3 g/s. Several 
differences are observed between the distributions. Roughly 50% of the predicted current distribution lies 
within the inlet channel, compared with 30% of the actual current. The predicted current distribution is 
predominantly radial, while the measured current contains a significant axial component, stretching from the 
cathode base toward the anode tip. Similar differences occur between measured and calculated current 
distributions in the cylindrical Stuttgart MPD thrusters 88 , and have been ascribed to variations in the cathode 
emission distributions between quasi -steady and steady-state operation 55,88 . Electrode processes have yet to be 
coupled with MPD flow models, however, and it is thus arguable whether the discrepancy between 
measurement and prediction lies in physical differences between quasi-steady and steady-state electrode emission 
processes, or in the inadequate modeling of relevant electrode phenomena. An experimental comparison of 
anode current distributions for a segmented anode thruster operated under quasi-steady and steady-state 
conditions would help to resolve this issue, and contribute to our understanding of whether quasi-steady MPD 
thruster data accurately reflects steady-state operating conditions. 

Preliminary modeling results have been reported on the influence of the electrode geometry 43 ** 7,52,56,95 on 
the self-field thruster current distribution. In general, the analyses indicate that contouring the electrode 
geometries may be an effective way to control the current distributions. Converging-diverging anode geometries 
appear to provide a more uniform current distribution along the electrodes than constant area channels, where 
the current tends to congregate at the inlet and exit regions of the thruster. The more uniform current 
attachment achieved with the converging-diverging anode geometries may help to mitigate the localized heating 
and associated erosion of the thruster electrodes. 

Voltaze. 


In addition to anode current fractions, Table 3 displays the measured and predicted discharge voltages 
for the University of Stuttgart ZT-1 thruster. The predicted voltage drops were calculated by integrating the 
radial electric field from cathode to anode, excluding electrode fall voltages or sheath potentials. The Stuttgart 
model predicts a potential drop across the plasma of 14 volts, which is in fair agreement with the total discharge 
voltage of 19 volts. The LeRC model, however, predicts a potential drop of only 3.5 volts across the plasma, 
much lower than either the Stuttgart model prediction or the measured total voltage. The calculated plasma 
potential is inversely proportional to the electrical conductivity, and the discrepancy in the model predictions 
arises from differences in the electrical conductivities calculated in each model. The contribution of electrode 
fall and sheath voltages to the total discharge voltage was not measured in the steady-state device, hence it is 
uncertain which model more accurately reflects the actual plasma potential. Experimental measurements of 
anode fall voltages or plasma potentials, in conjunction with anode current distribution measurements, would 
allow significant refinement of current MPD thruster models and would be of tremendous aid in the 
development of future electrode sheath and boundary flow models. 

Figure 13 displays measured total 87 and calculated plasma 62 voltage-current characteristics for the 
Princeton HSFAT for a variety of currents and mass flow rates. The model does not include electrode fall or 
sheath potentials, hence predicts a voltage substantially below the total measured discharge voltage. 

Nevertheless, the calculated plasma drop increases with increasing current, though at a slightly reduced rate 
compared with the total discharge voltage. In addition, the calculated potentials reproduce the dependence on 
mass flow rate observed in the total voltage measurements. 

Niewood and Martinez-Sanchez 63 have recently modeled the effects of high Hall currents near constant 
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area MPD thruster anodes. In addition to highly skewed current contours (Fig. 8), the Hall effect leads to 
anode starvation and a subsequently large anode potential drop. At a discharge current of 30 kA and a mass 
flow rate of 4 g/s, the calculated anode fall was approximately 19 volts, out of a total calculated discharge 
voltage of 26 volts. Heimerdinger 52 measured a near-anode fall of 25 volts for a fully flared thruster operating 
under similar conditions, with a total voltage drop of approximately 40-50 volts. While the model predictions 
are low compared to similar experimental measurements, they indicate that at least under the reported operating 
conditions Hall currents may have a significant impact on the anode fall voltages. Such preliminaiy results 
should be followed up with additional experimental measurements and modeling predictions under a variety of 
operating conditions to establish the impact of electrode-boundary Hall terms on self-field MPD thruster 
performance. 


Applied-fi e ld M cxfels, 

The added complexity of applied magnetic field terms in the MPD acceleration process and the 
abundance of self-field experimental measurements have kept the numerical modeling efforts focused primarily 
on self-field MPD thrusters. However, interest has been rekindled in applied-field modeling efforts by the 
resurgence of applied-field MPD thruster experiments. 

One of the earlier treatments incorporating applied magnetic field effects with MPD thruster flow 
properties was published by Krulle 42 , who assumed an applied-field distribution and used simplified flow 
equations to study electromagnetic acceleration mechanisms at low P/m values. It was assumed that the plasma 
was fully ionized, and the self-induced magnetic field was negligible in comparison with the applied magnetic 
field. For the operating regimes considered (Pi m < 1 x 10 J0 A 2 -s/kg), it was determined that a large fraction 
of the total thrust was due to a pressure force arising from a radially compressive j x B force acting on the 
plasma near the cathode tip. The model also showed that substantial current was carried out of the thruster and 
into the plume region as the applied magnetic field is increased, in qualitative agreement with experiment 89 . 

Fradkin 89 ’ 90 performed experimental measurements on an applied-field MPD thruster in the 10-35 kW 
operating range, and derived simple analytic expressions for the voltage, thrust, thrust efficiency, and thermal 
efficiency which were in reasonable agreement with the experimental results. Fradkin attempted to explain his 
measurements using the applied field torque model discussed in Nerheim and Kelly 9 . His model was 
incomplete, however, in that it did not include recovery of thermal energy during the expansion through the 
magnetic nozzle. 

Tanaka and Kimura 54 analyzed applied-field acceleration mechanisms using two-dimensional 
electromagnetic field equations and quasi- ID fluid equations with a number of additional simplifying 
assumptions. They found that current distributions along the anode were strongly influenced by variations in the 
applied-field strengths, while current distributions along the cathode shifted from the base toward the cathode tip 
as the applied-field strength increased. Significant Hall currents were produced by the interactions of the 
discharge current with applied magnetic field. The Hall currents interacted with the applied magnetic field to 
produce plasma rotation, which was converted into axial momentum as the plasma expanded through the 
magnetic nozzle. The addition of joule heating due to the applied-field acceleration mechanism was found to be 
negligible for the operating regimes considered (discharge currents of 1000-2000 A, propellant mass flow of 0. 1 
g/s argon, axial magnetic field strengths of 0. 1-0.2 T at the cathode tip). 

Plasma flow through magnetic fields has long been of interest to workers in the plasma fusion energy 
program, whose somewhat diametric goal is to magnetically confine the plasma for as long as possible. The 
expertise in plasma/field interactions gained by researchers in the magnetic fusion community has recently been 
applied to studies of high power, high thrust, applied-field plasma thrusters for spacecraft propulsion 91 . Gerwin 
et al. n have modeled the flow of an ideal MHD plasma through a magnetic nozzle for a wide range of plasma 
densities and temperatures, and provide a detailed discussion of issues related to ideal plasma acceleration and 
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detachment from magnetic field lines. Detachment from magnetic field lines in the absence of meridional 
electric currents is analyzed by Hooper 93 , who finds that separation is significantly constrained unless electric 
currents flow across the magnetic field. 

Plasma Instability Models. 

Models have been developed to investigate the causes of both macroscopic and microscopic plasma 
instabilities associated with MPD thruster operation. As noted in previous sections, the severe electrode erosion 
observed at high values of J 2 /m corresponds to the onset of voltage oscillations and unsteady thruster operation. 
Several mechanisms have been proposed to explain the onset of these oscillations, including anode mass 
starvation®, flow choking due to enhanced back-EMF* 9 *™, and the triggering of electrothermal and gradient- 
driven instabilities as the plasma approaches full ionization 13,71,72 . It appears that different operating conditions 
may trigger one or a combination of the proposed mechanisms, thus limiting thruster performance and lifetime. 

On a less dramatic scale, plasma microinstabilities may lead to increased plasma resistivity, leading to 
an increase in the voltage drop and decreasing the thruster efficiency. Rempferer et a/. 73 used a linear 
dispersion relation to show that electron acoustic wave instabilities are possible in MPD thrusters under normal 
operating conditions. Choueiri et a/. 20,74,73 derived general nonlinear dispersion relations for current-driven 
plasma instabilities, and determined that lower-hybrid waves are the dominant instability for most MPD thruster 
operating regimes. As noted in a previous section, the existence of these waves has recently been verified in 
both low power and high power MPD thruster experiments 20 ' 25 . 

Plasma instability models have so far been uncoupled from the plasma flow models discussed earlier. 
Niewoood et al 94 considered the effects of microinstabilities on the thruster operating conditions using a quasi- 
ID fluid model 56 and a nonlinear dispersion relation for the modified two-stream instability. Plasma parameters 
obtained from the fluid model are used to determine whether the instability exists in a given region; if so, the 
anamolous collision frequency determined from the dispersion relation is used to determine new values of the 
ion and electron heating rates, which in turn affect the plasma flow model. The results of this iterative 
modeling show a slightly increased potential difference across the electrodes, with little noticeable change in the 
overall thruster performance. It would be of considerable interest to expand this effort by coupling other 
plasma microinstability models (particularly the lower hybrid instability) with a more detailed plasma flow 
model. This iterative approach between microinstability theory and macroscopic fluid model could provide 
substantial insight into the effects of plasma instabilities on MPD thruster performance. 

Conclusions 

Significant progress has been made in the establishment of a broad performance and lifetime database, 
and in the development of numerical models for MPD thrusters. Thrusters have been operated at specific 
impulses up to 7000 sec, at efficiencies approaching 40% and at power levels over 0.5 MW. Efforts directed at 
establishing empirical scaling laws for applied-field thrusters, identifying primary causes of thruster inefficiency, 
and establishing reasons for observed performance limitations have started to produce results which can be used 
to directly improve thruster performance. The observed correlations between the anode power loss, the thruster 
geometry, and the applied-field strength indicate that significant improvements can be obtained using appropriate 
anode shapes and propellant injection. The identification of the lower hybrid instability as the dominant source 
of plasma microtorbolence in thrusters operated over three orders -of-magnitude in power level will not only aid 
in the development of appropriate transport models, but also has significant implications for self-field thruster 
scaling. 


For the first time in the last 20 years cathode erosion rates have been measured which, assuming they 
are constant, are low enough to achieve the required lifetimes for the missions of interest. Identification of 
surface temperature, propellant purity, and chamber pressure as the controlling factors for cathode erosion will 
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greatly aid in demonstrating long life thrusters. However, the inherent coupling between MPD thruster 
performance and lifetime, as exemplified by the interplay between thruster geometry, chamber pressure, and 
cathode material backscattering rates, has been inadequately explored. 


Numerical models are approaching the degree of completeness required to obtain accurate predictions. 
Self-field MPD thrust values have been successfully calculated for a variety of thruster geometries over a wide 
range of operating conditions. Calculated anode current distributions agree well with steady-state thruster 
measurements, although major discrepancies are observed between calculated current profiles and measured 
interelectrode current distributions in quasi-steady thrusters. This disagreement between predicted and measured 
values may result from differences in cathode emission processes between quasi-steady and steady -state thruster 
operation, or may reflect the omission of relevant electrode processes (sheaths, fall potentials) in present MPD 
thruster models. Voltage trends are in general agreement with experimentally observed voltage-current 
characteristics, although again the lack of an accurate treatment of electrode fall phenomena in the fluid models 
precludes an exact prediction of the total discharge current. Applied-field thruster flow models are less 
developed than their self-field counterparts, although the effects of magnetic nozzles on flowing plasmas are 
under active investigation. Models have been developed to evaluate plasma microinstability effects in MPD 
thrusters, and to predict conditions leading to the onset of thruster voltage oscillations and subsequent electrode 
erosion. While significant work remains to be done, a basic understanding of the " full potential and limitation" 
of the MPD thruster is progressing via the combined development of experimental and numerical techniques. 
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Table 1: Recent Measurement of MFD Cathode Erosion Rates 


Cathode 

Impregnant" 

Propellant 

Mode of 
Operation” 

Erosion Rate 
(Mg/C) 

Comments 

Ref. 

2% ThOj 

Ar 

s.s. 

(1.5-5.4)xl0‘ 3 

U. C. P. 

18,33 



s.s. 

1.8X10' 1 

Low power (< 30kW) 

17 



s.s. 

3.0xicr 2 


31 

1% Th02 


s.s. 

(1.6 - 4.9)xl0' 3 

U. C. P. 

35 

2% ThOj 


s.s. 

1.9 x 10- 3 

U. C. P. 

33 


N, 

s.s. 

(5.0-6.0)xl0-* 

U. C. P. 

18 


h 2 

s.s. 

(6-20)xl0“* 

U. C. P. 

18 


N 2 +2Hj 

s.s. 

1.3xlO' 3 

U. C. P. 

18 

? 


Q.S.S. 

6.0x10-' 

25 days accumulated life 

39 



Q.S.S. 

8.0- 13.0 

Applied B Field 

23 

BaO 

Nj+2H 2 

Q.S.S. 

(2-4)xl0-‘ 

Self and Applied B Field 

29 

2%Th0 2 

At 

Q.S.S. 

1-60 


36 


* Base material is tungsten 
~S.S. = steady state 
Q.S.S. — quasi-steady state 
U.C.P. = ultra clean propellant 
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Table 2. Recent MPD thruster models. 


ZT-1 THRUSTER 


U. STUTTGART 
EXPERIMENT I MODEL 


NASA LERC 
MODEL 


ANODE CURRENT FRACTION 


1 

0.46 

HU 

0.47 

2 

0.27 


0.23 

3 

0.27 

wESm 

0.30 


VOLTAGE (V) 


DISCHARGE VOLTAGE (V) 


19 V 


14 V f 


Table 3. Current distribution and discharge voltage comparisons for the ZT-1 MPD thruster operating at 
6000 A, 6 g/s. Calculated voltages do not include electrode fall potentials. 


25 






































c. Half-scale flared anode (HSFAT) quasi-steady 
thruster run at power levels up to 5 MW (Ref. 87). 


Figure 1. Schematics of typical MPD thrusters tested since 1987. 
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Figure 2. Peak efficiency - specific impulse values reported since 
1987. 



Power, kW 


Figure S. Peak efficiencies vs. thruster power level for data reported 
since 1987. 
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Figure 4. Peak specific impulse vs. thruster power level for 
measurements reported since 1987. 
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Figure 5. Comparison of performance levels obtained using hydrogen 
and argon propellants. Thruster geometry shown in Fig. la, propellant flow 
rate of 25 mg/s, discharge current of 750A. 
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Anode Power/Total Power 



Applied Magnetic Field, T 


Figure 6. Anode power fraction and total anode power vs. applied magnetic 
field strength. Thruster geometry shown in Fig. la, 0.1 g/s argon propellant, 
discharge current of 1000A. 



Figure 7. Erosion rate as a function of the ambient pressure of a 2 % 
thoriated tungsten cathode in a nitrogen atmosphere (Ref 35). 
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Figure 8. Theoretical current contour lines showing strong Hall effect 
near electrodes. Planar geometry thruster, 31 kA, 4 g/s argon propellant. 
Adapted from Ref. 63. 


halt scale flared anode thruster 

THRUST CHARACTERISTICS 



Figure 9. HSFAT thrust characteristics, argon propellant (Ref. 62). 
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Figure 10. HSFAT thrust coefficients, argon propellant (Ref. 62). 
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a. Enclosed current predictions from Ref. 88. 


b. Enclosed current predictions from Ref. 62. 


Figure 11. Predicted enclosed current contour comparisons for the ZT-1 thruster 
operating at 6000 A, 6 g/s, argon propellant. Contour levels every 10%. 
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a. Predictions from Ref. 62. 



Figure 12. Comparison of predicted steady-stated and measured quasi-steady current 
distributions. HSFAT thruster, S g/s argon, 17.8 kA discharge current. 
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VOLTAGE (V) 



Figure 13 . Measured total voltage and predicted plasma voltage vs. current for 
the HSFAT thruster. Argon propellant. 
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